Disorder induced cross-over effects at quantum critical points 
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Critical properties of quantum spin chains with varying degrees of disorder are studied at zero 
temperature by analytical and extensive density matrix renormalization methods. Generally the 
phase diagram is found to contain three phases. The weak disorder regime, where the critical be- 
havior is controlled by the fixed points of the pure system and the strong disorder regime, which is 
attracted by an infinite randomness fixed point, are separated by an intermediate disorder regime, 
where dynamical scaling is anisotropic and the static and dynamical exponents are disorder depen- 
dent. 
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Quantum phase transitions take place at zero temper- 
ature due to quantum fluctuations by varying a quan- 
tum control parameter which could be either the 
density of impurities (for the metal-insulator transition 
1^ in three dimensions (3d) or for the superconducting- 
insulator transition |^ in 2d) or an external thermody- 
namical parameter such as the strength of the magnetic 
field in quantum spin glasses Q or the pressure in BaVS^ 
l^. The effect of quenched, i.e. time-independent disor- 
der at quantum phase transitions is the subject of recent 
intensive research ||^; we mention, for instance, the re- 
cent experimental studies about doping effects in quan- 
tum spin chains j?] and ladders [|j and a theoretical ex- 
planation about the non-Fermi liquid behavior in U and 
Ce intermetallics 

From a theoretical point of view several analytical and 
numerical studies have been performed mainly in low- 
dimensional systems j^]. In particular in one dimen- 
sion many, presumably exact, results have been obtained 
by the application of a strong disorder renormalization 
group (SDRG) method, originally introduced by Ma, 
Dasgupta and Hu pO|. For the random transverse-fleld 
Ising model Fisher |0| has shown that at the critical 
point the probability distributions of the couplings and 
the transverse fields broaden without limits under renor- 
malization, and as a consequence the SDRG becomes 
asymptotically exact at the so-called infinite randomness 
fixed point (IRFP). 

In the IRFP dynamical scaling is strongly anisotropic, 
such that the relaxation time, tr, and correlation length, 
^, are related as: Intr ~ ^/^. Scaling of the static singular 
quantities can be conveniently described by the magne- 
tization scaling dimensions, x, (and similarly for surface 
spins by x^) and by the critical exponent of the average 
correlation length, i^, which are given by [0|: 



The asymptotically exact SDRG analysis of the random 
transverse-field Ising model has been extended into the 
off-critical region i.e. into the Griffiths phase jl^] 
and the dynamical exponent, z, connecting time- and 
length-scales as tr ~ C^i has been analytically calculated. 

The SDRG method has been applied for a series of ran- 
dom quantum systems at the transition point ]To| , p^ -pT[ . 
Despite the fact, that in the absence of disorder, these 
systems have different type of (spatial and temporal) cor- 
relations, thus they are in different universality classes, 
the SDRG analysis indicates that if the disorder is strong 
enough all models are attracted by the IRFP of the ran- 
dom transverse-field Ising model. In that respect quan- 
tum fluctuations seem to be irrelevant, so that solely the 
strong disorder effects dominate the critical behavior. 

However the SDRG approach does not provide any in- 
formation on the behavior of the system at weak disorder. 
In this limit one generally invokes the Harris criterion 
M,E9[ according to which weak disorder is irrelevant if: 
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where vq is the correlation length critical exponent of 
the pure system. If vq < 2/d, thus disorder is relevant 
the IRFP could be strongly attractive as for the random 
transverse- field Ising model, while ii Vq > 2/d only a 
finite amount of disorder may bring the system to the 
IRFP. 

The aim of this Letter is to investigate in detail the 
weak to strong disorder crossover behavior at a quan- 
tum critical point and to present a possible common sce- 
nario. In view of the generality of this phenomenon and 
its possible relevance in realistic quantum systems it is 
of importance to perform a detailed study of the singular 
behavior of random quantum systems in the whole range 
of strength of disorder. 

As a subject of our study we have chosen two one- 
dimensional quantum models, the q-state quantum clock 



1 



model (CM) and the quantum Ashkin- Teller (AT) model. 
Our selection is motivated by the following facts: i) both 
models are self-dual, thus the location of the critical point 
is exactly known, ii) The SDRG analysis can be per- 
formed for both models leading to IRFP behavior for 
strong disorder, iii) Critical properties of the pure mod- 
els are known in details, in particular one can locate the 
part of the phase diagram, where the Harris irrelevance 
criterion is valid in the form of Eq.(|^). 

We start to introduce the quantum CM which is 
defined by the Hamiltonian: 



CM 
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in terms of a discrete spin variable, Si = 1, 2, . . . , g, at 
a lattice site, I, and M;''^|s;) = js; ± l,mod q) is a spin 
raising (lowering) operator. 

Similarly the quantum AT model is defined by the 
Hamiltonian lElll: 
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in terms of two sets of Pauli matrices, cr^'^' ,Tf'^' . For 
both models the couplings, Ji, and the transverse fields, 
hi^ are independent random variables, while for the AT 
model the coupling between the two Ising models, e, is 
disorder independent. Both models are self-dual, which 
amounts to the invariance of the Hamiltonians in Eqs.(|^) 
and (^) under the transformation Ji ^ hi. 

The pure CM, with Ji — J and hi — h, ioi q > -i 
has an extended critical phase at both sides of the self- 
duality point, which is located at J — h. According to 
our numerical results, obtained by the density matrix 
renormalization group (DMRG) method |^ the critical 
exponents of the q = b model at the self-duality point 
are given by: 



xcM = 0.105(5), x'cM = 0.23(1). 



(5) 



Since the self-duality point is located in the middle of a 
critical phase vgm is formally infinity, thus small disorder 
is irrelevant, according to Eq.(||). 

The pure quantum AT model is critical along the self- 
duality line for — 1 < e < 1 with the critical exponents: 
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which are conjecturedly exact [^,^. For — 1/\/2 < e < 
— 1/2 the correlation length exponent of the pure model, 
i^AT, exceeds the value of 2, thus the Harris irrelevance 



criterion in Eq.(||) sets in. (We note that i^at stays for- 
mally infinite in the whole region of -1 < e < -l/\/2, 
which is the central part of the so called critical fan p]|.) 

The SDRG transformation can be performed for both 
models leading to similar decimation rules: 
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where the first (second) equation refers to the elimina- 
tion of a strong bond (field) of strength, Ji {hi). For 
the CM K = [1 - cos(27r/g)]/[l -I- 52,q] |l| and q does 
not renormalize under the transformation. On the other 
hand for the AT model, where, k = 1 + e, the coupling, 
e, enters into the renormalization as e = e^(l -I- e)/2. The 
RG equations in Eq.(^ are very similar to that of the 
random transverse-field Ising chain, which is recovered 
with K, — 1. Both systems will scale to IRFP of that 
model for strong enough initial disorder and for k > 0. 
Note, however, that for < k < 1, i.e. for g > 4 for the 
CM and — 1 < e < for the AT model, in some cases 
the generated new couplings/fields could be larger than 
the decimated ones, therefore the SDRG is not valid for 
weak disorder. 

To see in detail the weak-to-strong disorder cross-over 
phenomena we perform a parallel numerical study about 
the critical behavior of the two models with varying de- 
gree of disorder. We use the DMRG method and cal- 
culate the average magnetization profiles, to;, fixing the 
spins at one boundary, si = 1, and erf = = 1, respec- 
tively, whereas the spins at the other end of the chains, 
at I — L, are free. In particular we consider the scaling 
of the magnetization for spins at the bulk, when I = L/2 
and nib = '^1^/21 ^ind that at the surface, when iris = rriL: 
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from which the scaling dimensions, x and are calcu- 
lated. 

Useful information about the dynamical exponent, z, 
can be obtained from the distribution function of the sur- 
face magnetization, which has the limiting behavior : 



P(lnTOs) m^J^ , 



. 
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The scaling exponent, z', in Eq.(^ should be compared 
with that of due to pure quantum fluctuations, = 1, 
and then z — max(l, z'). 

For both models we perform the calculations at the 
self-duality point, which corresponds to the critical point 
of the pure systems, as well as, according to the SDRG 
analysis, to the fixed-points of the strongly disordered 
systems. By choosing the same distribution function, 
P{y)dy, of the couplings J; and that of the fields hi, 
we ensure to stay at the self-duality point and use the 
parametrization: 



P{y) = Ay 



-1+1/A 



0<y<l 0<A<oo. (10) 
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One can see here that A plays the role of the strength of 
disorder: the pure system is recovered as A ^ 0, whereas 
as A — > oo the distribution corresponds to that in the 

iRFP fnti. 

In the actual calculations we used q = 5 for the CM 
while for the AT model we took e = —0.75. The bulk 
and surface magnetizations are calculated on relatively 
large finite systems up to L = 32 and 10000 disorder 
realizations were used, except for the largest sizes of 
the AT model, when we have some 1000 samples. For 
stronger disorder, i.e. for larger A the energy gap in a 
typical sample is decreasing, therefore to maintain suffi- 
cient numerical accuracy of the code one needs to keep 
relatively more states in the DMRG procedure. As a 
consequence the computational demand of the calcula- 
tion is strongly increasing with A, therefore we have re- 
stricted ourselves to calculate only at few points in the 
strongly disordered regime. Generically the DMRG al- 
gorithm performs rather well for disordered systems, as 
it is known also from other several examples MjlJ]. 
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FIG. 1. Plots of 1/z' vs. A for the random AT model with 
e = -0.75 (■) and for the CM with <j = 5 (o). The limiting 
values of Aq and Aoo, where 1/z' = 1 and 1/z' — 0, respec- 
tively are model dependent. The error bars in the figure are 
due to finite size efi'ects and disorder averaging. 

Figure |l| shows the inverse of the disorder induced dy- 
namical exponent, 1/z', as calculated from Eq.(||) for dif- 
ferent values of the disorder strength, A. In each cases 
1/z' is found to be a monotonically decreasing function 
of A, such that 1/z' is lowered below the critical limit of 
1/z' = 1/zq = 1 for a model dependent disorder strength, 
Aq. Then, by further increase of A we obtain for the true 
dynamical exponent of the model, z = z' , to diverge, thus 
l/z = in the range of Aco < A < oo, where the upper 
disorder limit, Aoo is also model dependent. 

Next in Fig. ^ the scaling exponents of the average 
magnetization, x and X ^ clS calculated from Eq.(H) are 
shown for different strengths of disorder. Apparently, 



for weak enough disorder both the bulk and the surface 
exponents stay constant keeping the corresponding pure 
system value, which is in accordance with the Harris ir- 
relevance criterion. Then, by increasing the strength of 
disorder over a limiting value, which approximately cor- 
responds to Aq, found in Fig.l, the exponents start to in- 
crease. This variation is continuing until the exponents 
reach their IRFP values, as given in Eq.(|l|), then they 
stay constant for further increase of the disorder. This 
second limiting value of the strength of the disorder is 
approximately consistent with A^o, what can be identi- 
fied in Fig. |l] for the different models. Note, that the 
average behavior of the systems, as presented in Fig. 
is dominated by rare realizations, which are samples with 
mb{ms) = 0{1), but which occurs with a very low prob- 
ability. Therefore to obtain accurate averages one have 
to calculate on a large amount of samples, and the accu- 
racy of the calculation turned out to be larger for bulk 
exponents than that of the surface ones. 




FIG. 2. Scaling exponents of the average bulk and surface 
magnetization as a function of the strength of disorder for the 
q = 5 CM (o) and for the AT model with e = -0.75 (■). The 
limiting values of the exponents, corresponding to that of the 
IRFP are denoted by dashed lines. 

From the numerical and the SDRG results obtained 
for the two models we can extract a common route from 
the weak to strong disorder crossover, which is expected 
to be generally valid for a class of random critical quan- 
tum systems. This scenario is summarized in a generic 
RG phase diagram depicted in Fig. 3 in terms of the 
strengths of disorder. A, and a model dependent param- 
eter, say LJ, which is used to characterize points in a crit- 
ical line, (c.f. u! = e for the AT model, uj — A — q ioi the 
CM and ut — q/2 ~ 1 for the quantum Potts model). For 
Lu > the system is in the strong disorder (SD) phase, 
where the IRFP is strongly attractive. In the other part 
of the phase diagram, for a; < 0, there are two more 
phases: the weak disorder (WD) and the intermediate 
disorder (ID) regions. Here the SD part of the phase 
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diagram is attracted by the IRFP, while the critical be- 
havior in the WD regime is governed by the fixed points 
of the pure system, which are located at A = 0. In the 
ID regime, where there is a competition between quan- 
tum fluctuations and disorder effects, dynamical scaling 
is anisotropic, 1 < z < oo, and the static and dynamical 
critical exponents are disorder dependent. At the bound- 
aries of the ID region there are 1/z = and 1/z = 1, 
respectively, while at A = 0, i.e. in the pure system limit 
they are at w = and co — coq. In the latter case for the 
pure model with a parameter luq the Harris criterion in 
Eq.(||) is saturated, i.e. i'o(^o) — 2/(i. 
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FIG. 3. Schematic RG phase diagram of a quantum spin 
system having a critical line parametrized hy u, in the pres- 
ence of disorder of strength, A. For details see the text. 



We believe that the above scenario is relevant for some 
higher dimensional systems and therefore could have 
experimental relevance, too. As an example we men- 
tion the two dimensional random transverse-field Ising 
model where by a numerical implementation of the SDRG 
method an IRFP is found which is attractive for 
strong disorder both for spin-glasses (SG-s) and random 
ferromagnets (FM-s). In numerical studies, however, a fi- 
nite dynamical exponent, z, is obtained for the SG model 
p6[ , while for the random FM z is divergent at the critical 
point ||2^. A possible explanation of this controversy is 
that the studied SG is still in the ID region, whereas the 
random FM is in the SD phase. Further examples could 
be provided by random antiferromagnetic spin systems: 
chains, ladders and 2d models, where often conventional 
fixed-point behavior is obtained |l6yi7|j2^ . 
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